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Abstract 

An extremely accurate attitude determination 
has been developed by Aerojet ElectroSystems 
Company for space applications. This system 
uses a high sensitivity star sensor in which the 
photo multiplier tube is subject to noise generated 
by space radiations. The space radiation induced 
noise arises from trapped electrons, solar protons 
and other ionizing radiations, as well as from dim 
star background. The individual noise components 
have been successfully simulated for pre-flight 
evaluation of the star sensor performance. 

The steady trapped electron and the intermit- 
tent solar proton environments for a high altitude 
orbit are related to solar activity. The solar activ- 
ity and hence the electron and proton environments 
are predicted through the end of the twentieth cen- 
tury. The available data for the response of the 
phototube to proton, electron, gamma ray, and 
bremsstrahlung radiations are reviewed and new 
experimental data is presented. Prompt and de- 
layed pulses from the faceplate /photocathode and 
contributions from the dynodes are isolated, and 
pulse height distributions are given for each. 

The predicted electron and proton environ- 
ments are then combined with phototube response. 
The average noise due to the relatively steady 
electron environment is of the order of lO^pulses 
per second compared with 2 x 10^ pulses per 
second for the signal from a single 4.6 magnitude 
star. For a typical worst solar proton event, the 
noise will be 10^ to lO^pulses per second, which 
is comparable to the signal from a nominal 4.6 
magnitude star. This noise will typically last 
about one day. The model for the dim star back- 
ground averages 3 x 10^ pulses per second but 
varies widely. 

A simulation was developed which represents 
the characteristics of the effect of radiations on 
the star sensor, including the non- stationarity of 
the backgrounds. 

Introduction 

Aerojet ElectroSystems Company has develop- 
ed a high accuracy attitude determination system 
for space applications. The system contains a 


high sensitivity star sensor subject to space 
radiations. The object of this study is the pre- 
diction of the effects of those radiations. 

Star Sensor Description 

The optical portion of the star sensor collects 
light from stars and focuses it on the surface of a 
reticle. When a star image falls on a slit in the 
reticle, the light from that star can reach the 
photocathode of an EMR 541 E- 01 photomultiplier. 

The star sensor electronics were designed 
for photon counting. The counting method elim- 
inates the effects of stochastic variations in anode 
pulse amplitudes, i. e. , there is one count for 
each cathode photoelectron detected. The cir- 
cuitry which detects anode pulses outputs a nor- 
malized 50 nanosecond pulse for each anode pulse 
detected. These output pulses are counted 
throughout a 140 microsecond count cycle. The 
count value accumulated in each count cycle is 
transmitted to the ground for processing. The 
center of the star image requires about 230 
microseconds to traverse the slit, so the time 
at which the center of the image is at the center 
of a slit can be determined from the count values. 

When several photons reach the photocathode 
simultaneously, their anode pulses overlap, and 
there is only one normalized pulse for counting 
purposes. Thus when a high energy electron or 
proton causes the simultaneous emission of many 
protons, the star sensor electronics normally 
outputs only one normalized pulse. However, de- 
layed photon emissions from an incident electron 
or proton result in several output pulses. 

The counts from a star dimmer than seventh 
magnitude are comparable to the stochastic vari- 
ations in the total background signal and such 
stars cannot be detected. The average contribu- 
tion of stars dimmer than 4. 6 magnitude is ex- 
pected to be 4. 2 counts per count cycle although 
it is much higher in the galactic plane. 

Electron and Proton Environments 

The trapped electron and solar proton environ- 
ments are both capable of generating noise in star 
sensors. Both environments are affected by solar 
activity. 
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The detailed relation between solar eruptions 
and the solar proton and trapped electron environ- 
ments is very complex, and we are not currently 
able to predict satisfactorily the specific causitive 
solar eruptions. However, estimates of the ex- 
pected electron and proton fluxes can be related to 
expected smoothed sunspot numbers. 


F ( >E, >0 max )_ (2) 
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where 


The expected smoothed sunspot numbers are 
represented by 


R ~ R 0 + AR sin 2 



where 


R Q is the expected minimum R 
(typically 2-10) 


AR is the expected span of R 
(typically 50-200) 


( 1 ) 


R is the smoothed sunspot number 
E is proton energy (Mev) 

F is the expected number of events per year. 

For protons with energies around 50 Mev, solar 
proton events will generally last from 0. 5 to 
3 days, with one day being typical. Lower energy 
protons often arrive later and last for a longer 
total time. Duration expectancy may be scaled 
with >/e 7 Total proton fluence for all events in a 
given period of time may be expected to be 2 or 
3 times that from the expected single largest 
event during that period of time. 


t is the time from solar minimum 
(years) 

Equation 1 predicts a rise up to 6 months pre- 
maturely but is otherwise a good representation 
of expectancy. 

The dates of certain minima and predicted 
minima (Reference 1) are 


Beginning cycle 


18 

February 1944 

19 

April 1954 

20 

October 1964 

21 

December 1974 

22 

September 1985 

23 

January 1997 


Equation 2 represents a satisfactory correla- 
tion, based on available data from Cycles 19 and 
20. However, at energies less than 30 Mev, it 
tends to overpredict. It also tends to overpredict 
the number of smaller proton events. Recent data 
may be conveniently found in the issues of Refer- 
ence 2. 


For the trapped electron environment, the 
following model will be taken as representative 
of high altitude equatorial orbits: 



lxlO 7 e 


E 

0.21 


For the parameters in Equation 1, 


Cycle No. 

R o 

AR 

Estimated 6-year 
Mean of R Near 
Solar Maximum 

19 

6 

195 

150 

20 

8 

103 

85 

Average 

5 

100 

79 


Cycle 21 will likely be close to average. Cycles 
22 and 23 may be below average, though it is too 
early to predict with any confidence. For pre- 
dicting cycles beyond #20, the average values can 
be used. 


R/100 (3) 

E 

2 1 1 

+ 2x10 6 /_JL'\ e~ °* 50 electron 
\ 100 / cm 2 sec 

where 

R is the smoothed sunspot number 

E is electron energy (mev) 

A contribution to the electron flux below 0. 1 Mev 
has been neglected. 


For a high altitude orbit, solar proton events 
may typically occur during the order of 4% a year 
near solar maximum. Most of the time energetic 
solar protons will not be encountered. The ex- 
pected frequency of solar proton events larger 
than a given / max (proton/ cm 2 sec) may be 
taken as 


Equation 3 represents the mean rate averaged 
over one day. The diurnal extremes are typi- 
cally a factor of 2 higher at 0830 local time and 
a factor of 5 lower at 2400 local time. 

Primary cosmic ray fluxes will generally range 
from 4 to 9 particle/cm 2 sec, varying inversely 
with solar activity. They may be further de- 
pressed during solar proton events. The average 
energy is about 3 Bev. 
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Photomultiplier Irradiation Results 


The phenomenon of ionizing radiation inter- 
acting with a vacuum phototube device (e. g. , photo- 
multiplier) is believed to occur as follows: (1) a 

prompt pulse of 1 to several cathode photoelectrons 
equivalent followed by (2) a train of single photo- 
electron pulses with the rate of occurrence dying 
out with time (References 3 and 4). Pulse heights 
for the prompt pulse and also total pulses per in- 
teraction are shown in Table I, as deduced by 
various investigators for several vacuum photo- 
tubes and kinds of ionizing radiation.* 

The prompt pulse is believed caused by Ceren- 
kov radiation in the glass faceplate, while the de- 
layed pulse train is thought to be caused by excita- 
tions in the faceplate and subsequent decays with 
emission of optical quanta. Dynode pulses are re- 
latively few in proportion, being less than 0. 1% of 
the total for electrons, bremsstrahlung, and 
gammas and likely not greater than 0. 1% for pro- 
tons. For very large pulse sizes (above perhaps 
20 photoelectrons equivalent) the dynode pulses 
dominate. However, for photon counting star 
sensors, the dynode pulses can be neglected. 

Aerojet, in cooperation with others, has per- 
formed several experiments to measure the pulse 
height distributions from the photocathode and 
dynode s of photomultiplier tubes, due to electron, 
bremsstrahlung, and proton radiations. ** To ob- 
tain total radiation information, a photomultiplier 
tube may be used directly. The EMR 541E-01, 
with trialkali photocathode on 7056 glass (approx- 
imately S- 20 response) and copper -beryllium Ve- 
netian blind dynode s, is used in the star sensor. 
However, to separate photocathode and dynode ef- 
fects, it was necessary to find a device with large 
separation of photocathode from the dynode s. An 
ITT S-412 image dissector tube, trialkali photo - 
cathode on 7056 glass (approximately S-20 re- 
sponse) but with silver- magnesium box-and-grid 
dynode s, provided 4 inches separation and was 
selected for the supplemental experiments. The 
effective photocathode area of this tube is approx- 
imately 1.8 x 10 6 times smaller than for the EMR 
541 tubes, which further provides excellent isola- 
tion of the dynode effects. The faceplate is approx- 
imately 35% thicker. The wide separation of the 
aperture from the photocathode also permits di- 
rect determination of the pulse height distribution 
of delayed pulses without elaborate timing equip- 
ment, since prompt photocathode pulses are re- 
duced to single- electron events. 


From Sr-90/Y-90 beta and Co-60 and Ra-226 
gamma testing of the image dissector, two distinct 
dynode- region pulse height distributions were evi- 
dent. Whether the taller of two (Group A) origi- 
nates from the first dynode or the aperture plate 
was not determined. The delayed pulse height 
distribution for the photocathode region was de- 
termined from beta tests. Typical spectra for 
these tests are shown in Figure 1. 



Figure 1 90 

Pulae Height Spectrum Due To A Beta Source (Sr 


21 ) . 


Pulse height distributions from both photo- 
cathode and dynode regions were determined for 
protons at four energies from 22 to 39 Mev. A 
typical spectrum for these tests is shown in 
Figure 2. The response to protons is strongly de- 
pendent on proton energy, peaking very sharply 
around 34 Mev. 


* The luminescence parameter, L, is defined as the photocathode current divided by 
photocathode area, peak radiant sensitivity, and dose rate. 

L _ *c 

s U 0 ) f) ~c 

This luminescence parameter includes the Cerenkov radiation contribution. 

** The excellent work and cooperation of Dr. G.F. Knoll, University of Michigan; 

G. Bain and C.W, Freeborn, ITT; Dr. J. Ryan, MIT Lincoln Laboratory; and L.W. 

Morton, Aerojet, is hereby gratefully acknowledged. (Miss Freeborn is now at Aerojet.) 
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The photomultiplier was tested with Sr- 90 
bremsstrahlung and with bremsstrahlung from 
electrons at three energies from 0. 5 to 3 Mev to 
determine total pulses above a minimal threshold. 
The prompt pulse height distribution was not 
measured. This information has been inferred 
from data from other investigators and by analogy 
with the distribution for delayed pulses, 

Jerde, Peterson and Stein investigated the 
pulse height distribution for an RCA 7265 with 
S-20 response but a very thick soda-lime glass 
window. Their data for Co- 60 indicate about 2,5 
to 3 photoelectrons for the large- pulse height dis- 
tribution, which is to be associated with the 
prompt pulses, (Their reported value of 10 ap- 
pears applicable to only a few percent of the 
pulses in this distribution, and will therefore not 
concern us further here. ) Lacking specific data 
for the EMR 541E-01, we shall assume the mean 
prompt pulse size is 3 photoelectrons equivalent. 
Uncertainties due to difference in faceplate thick- 
ness make this estimate somewhere between good 
and up to a factor of 2 too high. 

Dressier and Spitzer offer measurements for 
the EMR 54 IE- 01 with Co-60 indicating approx- 
imately 20% of the total pulses are prompt. Jerde, 
et al. , report for the RCA 7265 about 3%, though 
their pulse height distribution data appear to indi- 
cate a somewhat larger value. Here we shall 
assume 20% for the EMR 541E-01. 


Jerde, et al. , for cosmic rays (predominately 
protons) report 10-20 photoelectrons per prompt 
pulse. This distribution was observed above 
5 photoelectrons pulse size and compares well 
with 13.8 photoelectrons per dynode pulse observed 
by AGC/ITT for the image dissector. Comparing 
with Figure 2, the cosmic ray pulses less than 3 
photoelectrons are coming predominately from the 
photocathode, while those greater than 5 photo- 
electrons are coming predominately from the dy- 
nodes. For large pulses, Cerenkov radiation 
would, therefore, appear to be a secondary effect 
with protons even up to cosmic ray energies. 

For single- electron counting. Aerojet observ- 
ed in the image dissector tube approximately 25% 
fewer pulses/rad with 39 Mev protons than with 
electrons and gammas. This is just about what 
would be expected if Cerenkov radiation is absent 
with protons, but all forms of ionizing radiation 
are equally effective in producing faceplate excita- 
tions. We hereby accept that this is true for 
photomultiplier tubes as well. Wolff (Reference 8) 
also finds that photocathode effects of image dis- 
sector tubes and photomultipliers are directly 
scalable with effective area. 

Zagorites and Lee (Reference 5) find radiation 
generated currents are essentially independent of 
gamma and X-ray energies except around 40-100 
kev, where it is probably enhanced by direct inter- 
actions with the photocathode or dynode s. Peak 
enhancement appears to be more than a factor of 3 
but not more than an order of magnitude. 

The rate of delayed pulses is reported by 
Dressier and Spitzer to decay in milliseconds. The 
effective time constant is believed to lengthen out 
somewhat in time due to the existence of a spec- 
trum of excited states. For most purposes, the 
early- time time constant will have the greatest 
impact. In order to have a number to work with 
tentatively use T = 1 x 10" 3 sec in the probability 
equation 

P = Ae t/r . (4) 

For most practical purposes, prompt pulses are 
instantaneous. 


Model for Photomultiplier Radiation Response 

A model pulse height distribution for the 
EMR 541E-01 irradiated by gamma rays, electrons 
and bremsstrahlung is shown in Figure 3. The 
models for delayed pulses and dynode pulses are 
based on experimental data, as described above. 
The assumed model for prompt pulses is, of 
course, dependent upon particle energy. For elec- 
tron energies below about 0. 5 Mev, gamma 
energies below about 0.8 Mev, and bremsstrahlung 
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Differential Fulsea per Unit Dose (Pulses/Rad-photoelectron) 


generated by electrons with energies below about 
1 Mev, the prompt pulses can probably be ne- 
glected* 



Figure 3 

tVlfe B&jghfc Spectra for An Em.541E-01 Photomultiplier (Co-60 Ganna Rays) 


The pulse height distributions of Figure 3 may 
be represented analytically: 


Photocathode — Delayed Pulses 


n ( h) = 8 x 10 ® exp /- 


h -1 




10. 873 J / 
rad-photoelectron, h < 3 


pulse/ 


= 4.2 x 10 12 x 10* h/0 ‘ 5 pul8e/ 
rad- photoelectron, h > 3 


(5) 


Photocathode — Prompt pulses 


n (h) ^ 0. 67 x 10 exp ( — 


h-3 


pulse/ 


2.62 

rad-photoelectron, h < 9 ( 6 ) 

~ 0. 35 x 10 12 X 10 _h/1 * 5 pulse/ 


rad-photoelectron, h > 9 


Integrated, the delayed pulses are about 

1 . 2 x 10 ^ pulse /rad, h > 0 . 1 ; 

1 . 0 x 10^ pulse/rad, h >0. 5; 
and the prompt pulses are about 

0 * 3 x 10^ pulse/rad, h >0. 1 or >0. 5 


) 4 exp — 

\ M 


Dynodes — Group A 
n(h) = 6 x 10 

photoelectron 

Dynodes — Group B 

n (h) = 1. 5 x 10~ exp 
photoelectron 

Integrated, the total pulses are 


-^pulse/rad- 


(7) 


H 


39/ 


Apulse/ rad- 


( 8 ) 


N a (>h ) = 8 x 10 4 x 10 ‘ h/3 pulse/rad (9) 
N 3 ( >h) = 6 x 10 5 x 10 " h ^°* ^ pulse/rad ( 10 ) 


For protons, the analytical representations are: 
Photocathode — Delayed Pulses 

N = 1,1 x 10’ pulse/rad (11) 

with distribution above h = 1 given by, 

( 12 ) 

N ( >h) = 1. 3 x 10 10 x 10" h/ 0 , 75 pulse/rad 


Prompt pulses likely occur only at relativistic 
proton velocities (i. e. , cosmic rays). 

Dynode — Group A (32 Mev protons) 


(■±-) 


n A I* 1 ) = 4. 4 x 1 0 3 exp 

rad-photoelectron 

Dynode — Group B (32 Mev protons) 


pulse/ 


n g (h) = 8 . 7 x 10 exp 


fzr) 


pulse/ 


(13) 


(14) 


rad-photoelectron 

Integrated these become 

N a (>h) = 6 x 10 4 x 10 pulse/rad (15) 

N b ( >h) = 4 x 10 5 x l(T h/ 10 ' 6 pulse/rad (16) 

Response to proton dose peaks sharply at 32 Mev 
and is about 4 Mev wide. 
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Alternatively, the integrated dynode responses 
may be represented by 


N a (. h) = 1.8 x 10 2 x io“ h/31 * 8 
pulse- cm 2 /proton 

Ng (>h) = 1.35 x 10~* x iQ-h/10.6 
pulse- cm 2 /proton 


(17) 


(18) 


In this case, virtually all noise generation will be 
by electron-gamma-delta ray (i. e. , Compton 
electron) interaction. Since only a small per- 
centage of incident electrons will yield delta rays 
exceeding the approximately 0. 5 Mev gammas, 
the mean prompt pulse height will be only h-1. 
Since many interactions will fail to produce 
Cerenkov radiation, it is possible that Equation ] 9 
overestimates the noise generation by as much as 
20 %. 


The above representations of dynode pulses 
are based on data for Ag Mg dynode s. Reed et al. 
data (Reference 6) indicate CuBe dynodes might be 
perhaps a factor of 3 less sensitive. Any such dif- 
ferences are overlooked here. 

The prompt pulses, including dynode pulses, 
should be assumed to occur in 10“ 7 seconds or 
less. The delayed pulse rate probability decays 
with a time constant of 10“ 3 seconds. For 
sapphire, Barr and Eberhardt (Reference 9) find 
the time constant becomes as long as several 
seconds to half a minute. 


Orbital Noise Prediction 


Therefore, the predicted mean orbital noise 
generation is 


N = 1.3 x 10 9 P 01 ^ 1 x 10~ 5 

rad sec 

= 1. 3 x 10^ pulse/ sec 


( 20 ) 


with virtually all pulses averaging h = 1.0. The 
expected number of pulses in one 140 millisecond 
counting interval is, therefore, 

N = 1.3 x 10 4 140 x 10~ 6 sec 

sec (21) 

= 1.8 pulse 


The electron environment will be present al- 
most continuously. Therefore, electron- induced 
noise will also be present almost continuously. 

The anticipated electron environment is given 
by Equation 3. One component is depressed by 
solar proton activity, while the other is enhanced. 
At R = 100, each component will yield about 10"° 
rad/ sec bremsstrahlung at the photomultiplier 
tube. The primary electrons below 5 Mev will all 
be stopped in the structure. The worst case for 
Cycles 20-21 will occur near solar minimum in 
1974- 1975. If it is assumed R ~10 is representa- 
tive of 1974, then noise in orbit will be generated 
by about 1 x 10~5 rad/ sec at the photomultiplier 

tube. 

If, from Table I, we accept 0.6 nanowatt 
(4100A) sec/rad cm for the luminescence of the 
EMR 541E-01-14 faceplate (h>0. 5) and 0.07 amp/ 
watt (4100A) peak radiant sensitivity, then the 
photomultiplier response will be 

N ( > 0.5) = 0. 6 x 10~ 9 ~ att sec 5cm 2 
rad cm 2 

(19) 

X 0.07 am P 1 P"! se 

watt 1.6 x 10’ 19 amp sec 

= 1. 3 x 10^ pulse/rad 


The high energy solar proton environment will be 
present only a few percent of the time. Over 90% 
of the year, solar proton- induced noise should be 
negligible. 

The expected solar proton environment is 
given by Equation 2. For a typical worst three 
years R-100. The star sensor structure, shield- 
ing, and surrounding hardware restrict protons 
reaching the photomultiplier to energies above 
50 Mev with an effective acceptance solid angle of 
about 0. 5 hemisphere = 0.25 sphere. For face- 
plate noise, the maximum solar proton environ- 
ment from Equation 2 is 

0 - 700 P roton / !2?\ 2 / 100 Mev \ 3 ^ 2 

max cm 2 sec ^100/ y 50 Mev J 

per year (3 years) (22) 

= 6000 proton/cm 2 sec (>50 Mev, 
omnidi r e ctional ) 
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Combining this environment with Equation 11, the 
faceplate noise will be 

N = 1. 1 X 10 9 6000 (0. 25) (1. 5) 

rad cm 6 sec 

V 1.8 X 10- 7 rad cm 2 (23) 

proton 

= 4. 4 x 10 5 pulse/sec 

where a spectrum correction factor F = 1. 5 has 
been assumed. This rate should endure for about 
one day. 

For the dynode pulses, it can be shown that 
the effective 32 Mev proton flux reaching the 
photomultiplier is about 

A(Z> max - 65 protons/cm 2 sec (24) 


Thus, using Equations 17 and 18, the two groups 
of dynode pulses will be 

. proton 


N, 


65 ^ (4 Mev) 

cm Mev sec 

2 pulse cm^ 

x 1.8 x 10"^ 10" 

proton 


(25) 

0.5/31.8 


= 1,2 pulse/ sec 


N b = 65 


proton 

7 

cm Mev sec 


(4 Mev) 


(26) 


i pulse cm^ .n c/io k 
x 1. 35 x 10- 1 1Q - 0.5/1 0.b 

proton 


= 8 pulse/ sec 


The dynode pulses represent only about 0. 03% of 
the total and are, therefore, negligible. 


Cosmic rays will generate not more than 


N = 1.5x109 g! 9 

cm* 1 sec 


(27) 


x 3 x 10' 8 p art i c i e - 4 pulse/sec 


rad cm^ 


including prompt pulses of perhaps h~7. Cosmic 
ray pulses also appear to be negligible compared 
to both the proton- generated and the electron 
generated noise. 

Thus, the maximum expected orbital noise is 
given by Equation 23, 

• 5 . 

N = 4. 4 x 10 pulse/ sec 


counting interval is therefore 


N = 4.4 x 10 5 PUlSeS 140 x 10 6 

SeC (28) 

= 61 pulses 

(For 1971-1973, these figures should be down by 
a factor of more than four. ) 

The pulse expectancies will not be entirely 
uncorrelated. (See Table I for the EMR 541E-01. ) 
After a faceplate interaction, whether producing 
a prompt pulse or not, the rate probability of de- 
layed pulses at time t afterwards is 


p (t) = = ~ 7 ~ exp ( * t/ ” > s < 29 > 

dt r 

where t may tentatively be taken to be 1 x 10“ J 
sec. Equation 29 applies to both electron- and 
proton- generated noise. Applied to a 140 micro- 
second counting interval, the probability of a 
pulse occuring is 

P(t) = 4(140 x 10-6 sec) ^ { _ t/10 _ 3 sec) 

1 x 10’^ sec 

- (30) 

= 0.56 exp ( -t/1 0~ sec) 

From Equations 21 and 28, it may be seen 
that the expectancy of prompt pulses in a given 
counting interval is 0. 36 and 12. 2 for the average 
electron and expected maximum proton environ- 
ments respectively. The expectancy of total 
pulses is 1.8 and 61 respectively. 


Star Sensor Simulation 


A computer program has been developed to 
simulate the star sensor performance. The out- 
put of the simulation program is used to provide 
inputs for development and qualification testing 
of the ground data processing. An earlier ver- 
sion of the simulation was used to evaluate the 
performance of the Kalman filter used in the atti- 
tude determination computations. The output 
from the star sensor simulation (star count pro- 
cessor) is a sequence of star sensor count values. 
Since the count cycle duration is 140 microseconds, 
almost one million of these values must be 
generated for every two minutes of star sensor 
operation simulated. 

The Star Count Processor is part of a com- 
prehensive simulation program modeling the en- 
tire satellite. Another processor in the simula- 
tion program is the Slit Crossing Processor 
which determines the times at which star images 
cross the center of a slit in the reticle. These 


with virtually all pulses averaging h = 1.0, The 
expected number of pulses in one 140 microsecond 
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slit crossing times are dependent upon the reticle 
slit geometry, the satellite dynamics being simu- 
lated, and star positions derived from a Star 
Catalog Data Set containing information from 
Reference 10, The output from the Slit Crossing 
Processor is a data set of slit crossing times 
along with the S. A. O. catalog number of the star 
and its visual magnitude. 

The Star Count Processor is capable of simu- 
lating a wide variety of background conditions 
specified by means of the input parameters dis- 
cussed below. For example, the parameter EP 
is used as the Poisson distribution parameter for 
modelling the arrival of prompt pulses caused by 
bremsstrahlung from trapped electrons. The 
parameter, EP, is the expected number of count 
cycles (140 microsecond intervals) between arri- 
vals of these prompt pulses. Since a Poisson 
process has exponentially distributed inter- 
arrival times (Reference 11), the interval between 
consecutive prompt pulses due to trapped electrons 
can be expressed as follows in terms of count 
cycles, 

DTNE = - EP In. (RU) (31) 

Where RU is a pseudorandom number gener- 
ated by the computer based upon a distribution 
uniform on the interval (0, 1). If we let STNE’ 
be the arrival "time" of the preceding prompt 
pulse due to bremsstrahlung, then the next prompt 
pulse arrives at 

STNE = STNE 1 + DTNE 

= STNE* - EP/n (RU) (32) 


Then attaching a prime to designate values of RE 
and RP from the preceding count cycle, we have 


RE = (RE 1 

+ n e EN ) DE 

(34) 

RP = (RP' 

+ rip PN ) DP 

(35) 


where 

n e = number of prompt pulses due to 

trapped electrons in the preceding 
count cycle 

rip = number of prompt pulses due to 
solar protons in the preceding 
count cycle 

When the processing for the current count 
cycle begins, its index, SNC, is obtained by add- 
ing one to the previous value of the index. If the 
previously generated value of STNE is less than 
or equal to the new value of SNC, then EN is 
added to RE*. A new value of DTNE is generated 
if STNE ^ SNC based on a new random number 
RU. The value of DTNE is then added to the cur- 
rent value of STNE. The process is repeated un- 
til the new value of STNE exceeds the current 
value for SNC. Similar logic is used for STNP 
with PN being added to RP. 

The star background effects vary as the star 
sensor optical axis sweeps across the celestial 
sphere. Let b Q be twice the number of radians 
swept in 1020 count cycles and K be the number 
of blocks of 1020 count values which have been 
generated so far in this simulation. The 
expression 


A similar logic applies to prompt pulses due to 
protons, either solar or cosmic, in terms of the 
input parameter SP, the expected number of 
count cycles between prompt proton-induced 
pulses. Hence 

STNP = STNP' - SP in (RU) (33) 

As noted previously, the delayed pulses fol- 
lowing a prompt pulse have a rate of occurrence 
dying out with time. This characteristic is 
modeled by means of the decay factors, DE and 
DP, for electrons and protons respectively. 

Then defining 

RE = expected number of output pulses 
in any count cycle due to electrons 

RP = expected number of output pulses 
in any count cycle due to protons 

EN = an input parameter related to ef- 
fect of each trapped electron 

PN = an input parameter related to ef- 
fect of each solar proton 


XS = 


.25 

1. 25 + sin b Q K 


(36) 


has a maximum value of 1. 0 when the sine func- 
tion is -1 and a minimum value of 1/9 when the 
sine function is +1. Then the expected contribu- 
tion from stars dimmer than seventh magnitude 
is 


BD = BDI * XS 


(37) 


where 

BDI = an input parameter equal to the ex- 
pected contribution of stars dimmer 
than seventh magnitude in the galac- 
tic plane. 


Contributions of stars brighter than 
seventh magnitude must be represented on an 
individual basis. Although stars dimmer 
than. 4. 6 magnitude could have been represented 
on the same basis as the brighter stars, it was 
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found that the extra computational cost was not jus- 
tified. Instead, the following simplified representa- 
tion was used. These intermediate intensity stars 
were simulated in terms of an exponential distribu- 
tion of interarrival times using the quantity SD as 
the expected interarrival parameter and 

SD = SDI/ (0.15 + XS) (38) 

where 

SDI = an input parameter 1.15 times the ex- 
pected interarrival time in count 
cycles in the galactic plane. 

The quantity SD is used in the same manner as the 
parameters EP and SP resulting in the equation 

STND = STND' - SD in(RU) + 3 (39) 

The three is added to avoid overlap of these stars. 

The intensity of the simulated star is given by 

QD = QM/ (.1 + KU) (40) 

where QM = input parameter nominally equal to one- 
tenth the response of the star sensor to a 4.6 magni- 
tude star. 

The timing of the slit crossing with respect to 
the count cycles is given in terms of another random 
number, Z, as follows 


RDl = (ZH-Z)QD 

(41) 

RD2 = QD 

(42) 

RD3 = (Z + ZH - 1.0 )QD 

(43) 


where ZH = an input parameter equal to the half 


3. In the single- electron counting mode, virtually 
all photomultiplier pulses originate in the face- 
plate/photocathode. This becomes important 
in design shielding. 

4. Star sensor response to the radiation environ- 
ment of space may be readily simulated with 
an appropriate computer program. 
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